Nondegenerate 3D complex Euclidean 
superintegrable systems and algebraic varieties 

E. G. Kalnins 
Department of Mathematics, University of Waikato, 
Hamilton, New Zealand. 
J. M. Kress 

School of Mathematics, The University of New South Wales, 
Sydney NSW 2052, Australia 
j.kress@unsw.edu.au 
and W. Miller, Jr. 
School of Mathematics, University of Minnesota, 
Minneapolis, Minnesota, 55455, U.S.A. 
miller@ima. umn.edu 

February 1, 2008 



Abstract 

A classical (or quantum) second order superintegrable system is 
an integrable n-dimensional Hamiltonian system with potential that 
admits 2n — 1 functionally independent second order constants of the 
motion polynomial in the momenta, the maximum possible. Such 
systems have remarkable properties: multi-integrability and multi- 
separability, an algebra of higher order symmetries whose representa- 
tion theory yields spectral information about the Schrodinger opera- 
tor, deep connections with special functions and with QES systems. 
Here we announce a complete classification of nondegenerate (i.e., 4- 
parameter) potentials for complex Euclidean 3-space. We characterize 
the possible superintegrable systems as points on an algebraic variety 
in 10 variables subject to six quadratic polynomial constraints. The 
Euclidean group acts on the variety such that two points determine 
the same superintegrable system if and only if they lie on the same 
leaf of the foliation. There are exactly 10 nondegenerate potentials. 
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1 Introduction 



For any complex 3D conformally flat manifold we can always find local co- 
ordinates x, y, z such that the classical Hamiltonian takes the form 

H = Tf— — \(Pi +P2 +pD + V ( X ,V, Z ), (x, y, z) = (xi,x 2 ,x 3 ), (1) 

i.e., the complex metric is ds 2 = X(x,y, z)(dx 2 + dy 2 + dz 2 ). This system 
is superintegrable for some potential V if it admits 5 functionally inde- 
pendent constants of the motion (the maximum number possible) that are 
polynomials in the momenta pj. It is second order superintegrable if 
the constants of the motion are quadratic, i.e., of the form 

S = J2 v)PjPi + W(x, y, z). (2) 
That is, {H, S} = where 

n 

{f,g} = ^2( d xJd P] g - d Pj fd Xj g) 

is the Poisson bracket for functions /(x, p), g(x, p) on phase space [U EJ [3j 
HI EJ El [3, [8] . There is a similar definition of second order superintegrability 
for quantum systems with formally self-adjoint Schrodinger and symmetry 
operators whose classical analogs are those given above, and these systems 
correspond one-to-one, [9]. (In particular, the terms in the Hamiltonian 
that are quadratic in the momenta are replaced by the Laplace-Beltrami 
operator on the manifold, and Poisson brackets are replaced by operator 
commutators in the quantum case.) Historically the most important su- 
perintegrable system is the Euclidean space Kepler-Coulomb problem where 
V = a/ \Jx 2 + y 2 -\- z 2 . (Recall that this system not only has angular mo- 
mentum and energy as constants of the motion but a Laplace vector that is 
conserved.) Second order superintegrable systems have remarkable proper- 
ties. In particular, every trajectory of a solution of the Hamilton equations 
for such a system in 6-dimensional phase space lies on the intersection of 5 
independent constant of the motion hypersurfaces in that space, so that the 
trajectory can be obtained by algebraic methods alone, with no need to solve 
Hamilton's equations directly. Other common properties include multisepa- 
rability (which implies multiintegrability, i.e., integrability in distinct ways) 
[H El El 13 El El El El EH HH H2] and the existence of a quadratic algebra of 
symmetries that closes at order 6. The quadratic algebra in the quantum 
case gives information relating the spectra of the constants of the motion, 
including the Schrodinger operator. 



2 



Many examples of 3D superintegrable systems are known, although, in 
distinction to the 2D case, they have not been classified, [131 HH EH EE [HI 
[18]. Here, we employ theoretical methods based on integrability conditions 
to obtain a complete classification of Euclidean systems with nondegenerate 
potentials. To make it clear how these systems relate to general second order 
superintegrable systems we introduce some terminology. A set of 2nd order 
symmetries for a classical superintegrable system is either linearly indepen- 
dent (LI) or linearly dependent (LD). LI sets can functionally independent 
(FI) in the 6-dimensional phase space in two ways: they are strongly func- 
tionally independent (FI-S) if they are functionally independent even when 
the potential is set equal to zero. They are weakly functionally indepen- 
dent (FI-W) if the functional independence holds only when the potential is 
turned on (example: the isotropic oscillator). Otherwise they are functionally 
dependent (FD). An LI set can be functionally linearly independent (FLD) 
if it is linearly dependent at each regular point, but the linear dependence 
varies with the point. An LI set can be FLD in two ways. It is weakly func- 
tional linear dependent (FLD-W) if the functional linear dependence holds 
only with the potential turned off and strongly functional linear dependent 
(FLD-S) if the functional linear dependence holds even with the potential 
turned on. Otherwise the set is functionally linearly independent (FLI). The 
Calogero and Generalized Calogero potentials are FD, and FLD-S, [9J. One 
property of FLD systems is that their potentials satisfy a first order linear 
partial differential equation. Thus they can be expressed in terms of a func- 
tion of only two variables. In that sense they are degenerate. This paper 
is concerned with a classification of functionally linearly independent poten- 
tials. As shown in [19], if a 3D second order superintegrable system is FLI 
then the potential V is must satisfy a system of coupled PDEs of the form 

V 22 = V u + A 22 V X + B 22 V 2 + C 22 V 3 , V 33 = V n + A 3 ^ + B 33 V 2 + C 33 V 3 , (3) 

V 12 = A 12 Vi + B 12 V 2 + C 12 V 3 , V 13 = A 13 Vi + B 13 V 2 + C l3 V 3} 

V 23 = A 23 Vx + B 23 V 2 + C 23 V 3 . (4) 

The analytic functions A % \ B % \ C l i are determined uniquely from the Bertrand- 
Darboux equations for the 5 constants of the motion and are analytic except 
for a finite number of poles. If the integrability conditions for these equa- 
tions are satisfied identically then the potential is said to be nondegener- 
ate. A nondegenerate potential (which is actually a vector space of potential 
functions) is characterized by the following property. At any regular point 
x o = (^0,2/0,^0), i-e., a point where the A^,B^,C l:i are defined and ana- 
lytic and the constants of the motion are functionally independent, we can 
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prescribe the values of V(x ), Vi(x ),V2(x ),V3(x ),Vii(x ) arbitrarily and 
obtain a unique solution of (jl]). Here, V\ = dV/dx, V 2 = dV/dy, etc. The 
4 parameters for a nondegenerate potential (in addition to the usual addi- 
tive constant) are the maximum number of parameters that can appear in a 
superintegrable system. A FLI superintegrable system is degenerate if the 
potential function satisfies additional restrictions in addition to equations 
(jlj). These restrictions can arise in two ways, either as additional equations 
arising directly from the Bertrand-Darboux equations or as restrictions that 
occur because the integrability conditions for equations (TjO) are not satisfied 
identically. In any case, the number of free parameters for a degenerate po- 
tential is strictly fewer than 4. In this sense the nondegenerate potentials 
are those of maximal symmetry, though the symmetry is not meant in the 
traditional Lie group or Lie algebra sense. Nondegenerate potentials admit 
no nontrivial Killing vectors. Our concern in this paper is the classification of 
all 3D FLI nondegenerate potentials in complex Euclidean space. In [20] we 
have begun the study of fine structure for second order 3D superintegrable 
systems, i.e., the structure and classification theory of systems with various 
types of degenerate potentials. 

Our plan of attack is as follows. First we give a brief review of the funda- 
mental equations that characterize second order FLI systems with nondegen- 
erate potential in a 3D conformally flat space. Then we review the structure 
theory that has been worked out for these systems, including multisepara- 
bility and the existence of a quadratic algebra. We will recall the fact that 
all such systems are equivalent via a Stackel transform to a superintegrable 
system on complex Euclidean 3-space or on the complex 3-sphere. Thus a 
classification theory must focus on these two spaces. Due to the multisep- 
arability of these systems we can use separation of variables theory to help 
attack the classification problem. In [2TJ we showed that associated with each 
of the 7 Jacobi elliptic coordinate generically separable systems for complex 
Euclidean space there was a unique superintegrable system with a separable 
eigenbasis in these coordinates. Thus the only remaining systems were those 
that separated in nongeneric orthogonal coordinates alone, e.g., Cartesian 
coordinates, spherical coordinates, etc. The possible nongeneric separable 
coordinates are known [22J so, in principle, the classification problem could 
be solved. Unfortunately, that still left so many specific coordinate systems 
to check that classification was a practical impossibility. Here we present 
a new attack on the problem, based on characterizing the possible super- 
integrable systems with nondegenerate potentials as points on an algebraic 
variety. Specifically, we determine a variety in 10 variables subject to six 
quadratic polynomial constraints. Each point on the variety corresponds to 
a superintegrable system. The Euclidean group E(3, C ) acts on the variety 
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such that two points determine the same superintegrable system if and only 
if they lie on the same leaf of the foliation. The differential equations de- 
scribing the spacial evolution of the system are just those induced by the Lie 
algebra of the subgroup of Euclidean translations. A further simplification 
is achieved by writing the algebraic and differential equations in an explicit 
form so that they transform irreducibly according to representations of the 
rotation subgroup SO (3, C ). At this point the equations are simple enough 
to check directly which superintegrable systems arise that permit separation 
in a given coordinate system. We show that in addition to the 7 superinte- 
grable systems corresponding to separation in one of the generic separable 
coordinates, there are exactly 3 superintegrable systems that separate only 
in nongeneric coordinates. Furthermore, for every system of orthogonal sepa- 
rable coordinates in complex Euclidean space there corresponds at least one 
nondegenerate superintegrable system that separates in these coordinates. 
The method of proof of these results should generalize to higher dimensions. 



2 Conformally flat spaces in three dimensions 

Here we review some basic results about 3D second order superintegrable sys- 
tems in conformally flat spaces. For each such space there always exists a lo- 
cal coordinate system x, y, z and a nonzero function X(x, y, z) = exp G(x, y, z) 
such that the Hamiltonian is ([T]). A quadratic constant of the motion (or gen- 
eralized symmetry) ^ must satisfy {H, S} = 0. i.e., 

-da 11 - G 2 a 2i - G 3 a 3i 
-da 1 ? - G 2 a 2 i - G 3 a 3 \ i ^ j (5) 
0, i, j, k distinct 

and 

3 

W k = \Y,a sk V s , fc = l,2,3. (6) 

s=l 

(Here a subscript j denotes differentiation with respect to Xj.) The require- 
ment that d Xe Wj = d x Wg, £ ^ j leads from to the second order Bertrand- 
Darboux partial differential equations for the potential. 

3 

[V SJ \a si - V se Xa SJ + V s {(Xa^ - (Xa sj ) e )} = 0. (7) 

s=l 

For second order superintegrabilty in 3D there must be five functionally 
independent constants of the motion (including the Hamiltonian itself). Thus 



■la 1 ; 1 + af 
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the Hamilton- Jacobi equation admits four additional constants of the motion: 

3 

S h = a lh)PkPj + W (h) = L h + W(h), h = 1, • • • , 4. 
j,k=i 

We assume that the four functions Sh together with H are functionally lin- 
early independent in the six-dimensional phase space. In [19j it is shown that 
the matrix of the 15 B-D equations for the potential has rank at least 5, hence 
we can solve for the second derivatives of the potential in the form If the 
matrix has rank > 5 then there will be additional conditions on the potential 
and it will depend on fewer parameters. D}J/\ + -D^V-j + D 3 s ^V 3 = 0. Here 
the A % \ B^ , C % i , are functions of x, symmetric in the superscripts, that 
can be calculated explicitly. Suppose now that the superintegrable system is 
such that the rank is exactly 5 so that the relations are only ([3]). Further, 
suppose the integrability conditions for system ([3]) are satisfied identically. 
In this case the potential is nondegenerate. Thus, at any point xo, where the 
A l i , , C u are defined and analytic, there is a unique solution V(x) with 
arbitrarily prescribed values of V"i(x ), V 2 (x ), V^xq), Vn(x ) (as well as the 
value of ^(xq) itself.) The points x are called regular. 

Assuming that V is nondegenerate, we substitute the requirement ([3]) 
into the B-D equations (I7|) and obtain three equations for the derivatives 
al . Then we can equate coefficients of Vi,V 2 ,V 3 ,Vn on each side of the 
conditions diV 2 3 = c^V 7 ^ = <9 3 Vi 2 , c^V^ = <9 2 V3 3 , etc., to obtain integrability 
conditions, the simplest of which are 

A 23 = B 13 = C 12 , B 12 -A 22 = C 13 -A 33 , B 23 = A 31 + C 22 , C 23 = A 12 + B 33 . 

(8) 

It follows that the 15 unknowns can be expressed linearly in terms of the 10 
functions 

A i2 , A 13 , A 22 , A 23 , A 33 , B 12 , B 22 , B 23 , B 33 , C 33 . (9) 

In general, the integrability conditions satisfied by the potential equations 
take the following form. We introduce the vector w = (Vi, V2, V3, Vn) T , and 
the matrices A^, j = 1, 2, 3, such that 

d Xj w = A®w j = 1,2,3. (10) 

The integrability conditions for this system are 

A? - Af = A® A® - A® A® = [A®, A®]. (11) 

The integrability conditions © and (ITT]) are analytic expressions in Xi,x 2 , x 3 
and must hold identically. Then the system has a solution V depending on 
4 parameters (plus an arbitrary additive parameter). 
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Using the nondegenerate potential condition and the B-D equations we 
can solve for all of the first partial derivatives a{ k of a quadratic symmetry to 
obtain the 18 basic symmetry equations, (27) in |19| . plus the linear relations 
(jHJ). Using the linear relations we can express C 12 , C 13 , C 22 , C 23 and -B 13 in 
terms of the remaining 10 functions. Each is a linear combination of the 
a with coefficients that are linear in the 10 variables and in the G s . 

Since this system of first order partial differential equations is involutive, 
the general solution for the 6 functions a jfc can depend on at most 6 param- 
eters, the values a jA; (xo) at a fixed regular point xo. For the integrability 
conditions we define the vector-valued function 

h(x,y,z) = (a 11 , a 12 , a 13 , a 22 , a 23 , a 33 ) T 

and directly compute the 6x6 matrix functions AP> to get the first-order 
system d Xj h = A^h, j = 1, 2, 3. The integrability conditions for this system 
are are 

A^h - Af ] h = A®A®h - A®A®h = [A®, A u) ]h. (12) 

By assumption we have 5 functionally linearly independent symmetries, 
so at each regular point the solutions sweep out a 5 dimensional subspace 
of the 6 dimensional space of symmetric matrices. However, from the condi- 
tions derived above there seems to be no obstruction to construction of a 6 
dimensional space of solutions. Indeed in [19] we show that this construction 
can always be carried out. 

Theorem 1 (5 ==>- 6) Let V be a nondegenerate potential corresponding to a 
conformally flat space in 3 dimensions that is superintegrable, i.e., suppose V 
satisfies the equations $31) whose integrability conditions hold identically, and 
there are 5 functionally independent constants of the motion. Then the space 
of second order symmetries for the Hamiltonian H = (p^+Py+pl)/ \(x, y, z) + 
V(x,y,z) (excluding multiplication by a constant) is of dimension D = 6. 

Thus, at any regular point (xq, yo, Zq), and given constants = a jfc , 
there is exactly one symmetry S (up to an additive constant) such that 
a h i(x , yo, z ) = a k K Given a set of 5 functionally independent 2nd order 
symmetries £ = {Sg : £ = 1, ■ • • 5} associated with the potential, there is 
always a 6th second order symmetry Sq that is functionally dependent on £, 
but linearly independent. 

Since the solution space of the symmetry equations is of dimension D = 6, 
it follows that the integability conditions for these equations must be satisfied 
identically in the a u As part of the analysis in reference [TH] we used the 
integrability conditions for these equations and for the potential to derive 
the following: 
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1. An expression for each of the first partial derivatives d^A^, diB^, c^C lJ , 
for the 10 independent functions as homogeneous polynomials of order 
at most two in the A i>j ' , C Vi ' . There are 30 = 3 x 10 such 
expressions in all. (In the case G = the full set of conditions can be 
written in the convenient form ( |59l) . ( 16T1) .) 

2. Exactly 5 quadratic identities for the 10 independent functions, see (31) 
in [19] . In Euclidean space these identities take the form 1^ — I^ e > in 
(|24|) of the present paper. 

In references [Tj5] we studied the structure of the spaces of third, fourth 
and sixth order symmetries (or constants of the motion) of H. Here the 
order refers to the highest order terms in the momenta. We established the 
following results. 

Theorem 2 Let V be a superintegrable nondegenerate potential on a con- 
formally flat space. Then the space of third order constants of the motion 
is 4- dimensional and is spanned by Poisson brackets Rjk = {Sj,Sk} of the 
second order constants of the motion. The dimension of the space of fourth 
order symmetries is 21 and is spanned by second order polynomials in the 6 
basis symmetries Sh- (In particular, the Poisson brackets {Rjk,S(} can be 
expressed as second order polynomials in the basis symmetries.) The dimen- 
sion of the space of sixth order symmetries is 56 and is spanned by third order 
polynomials in the 6 basis symmetries Sh- (In particular the products RjkReh 
can be expressed by third order polynomials in the 6 basis symmetries.) 

There is a similar result for fifth order constants of the motion, but it follows 
directly from the Jacobi identity for the Poisson bracket. This establishes 
the quadratic algebra structure of the space of constants of the motion: it is 
closed under the Poisson bracket action. 

From the general theory of variable separation for Hamilton- Jacobi equa- 
tions [221 [23] and the structure theory for Poisson brackets of second order 
constants of the motion, we established the following result [21]. 

Theorem 3 A superintegrable system with nondegenerate potential in a 3D 
conformally flat space is multiseparable. That is, the Hamilton - Jacobi equa- 
tion for the system can be solved by additive separation of variables in more 
than one orthogonal coordinate system. 

The corresponding Schrodinger eigenvalue equation for the quantum systems 
can be solved by multiplicative separation of variables in the same coordinate 
systems. 
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Finally, in [2T] we studied the Stackel transform for 3D systems, an in- 
vertible transform that maps a nondegenerate superintegrable system on one 
conformally flat manifold to a nondegenerate superintegrable system on an- 
other manifold. Our principal result was 

Theorem 4 Every superintegrable system with nondegenerate potential on 
a 3D conformally flat space is equivalent under the Stackel transform to a 
superintegrable system on either 3D flat space or the 3-sphere. 

3 Generic separable coordinates for Euclidean 
spaces 

Now we turn to the classification of second order nondegenerate superinte- 
grable systems in 3D complex Euclidean space. A subclass of these systems 
can be obtained rather easily from separation of variables theory. To make 
this clear we recall some facts about generic elliptical coordinates in complex 
Euclidean n space and their relationship to superintegrable systems with 
nondegenerate potentials (see [21] for more details). 

Consider a second order superintegrable system of the form H = J2k=i Pk+ 
V(x) in Euclidean n space, expressed in Cartesian coordinates x^. In analogy 
with the 3D theory, the potential is nondegenerate if it satisfies a system of 
equations of the form 



where all of the integrability conditions for this system of partial differential 
equations are identically satisfied, [251 HH]- There is an important subclass 
of such nondegenerate superintegrable systems that can be constructed for 
all n > 2, based on their relationship to variable separation in generic Jacobi 
elliptic coordinates. The prototype superintegrable system which is nonde- 
generate in n dimensional flat space has the Hamiltonian 



n 



V rj -V ll = Y J A m i^)V ll j = 2, ■ ■ ■ ,n, 



(13) 



n 



V kj = J2 Akj '*(x)Ve, l<k<j<n, 




(14) 
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This system is superintegrable with nondegenerate potential and a basis of 
n{n + l)/2 second order symmetry operators given by 

2 2 A 2 x<1 ' x<2 

Pi Pi ~\~ OlX^ + 2> Jij \%iPj •EjPi) Pi 2 Pj 2' ^ 7^ 3' 

Although there appear to be "too many" symmetries, all are functionally 
dependent on a subset of In — 1 functionally independent symmetries. A 
crucial observation is that the corresponding Hamilton- Jacobi equation H = 
E admits additive separation in n generic elliptical coordinates. 

x\ = c 2 WJ =1 (uj - ei)/Tl k ^i(e k - e*) 

simultaneously for all values of the parameters with e« 7^ tj if i 7^ j and — 
1, • • • , n. (Similarly the quantum problem if \l/ = is superintergrable and 
admits multiplicative separation.) Thus the equation is multiseparable and 
separates in a continuum of elliptic coordinate systems (and in many others 
besides). The n involutive symmetries characterizing a fixed elliptic separable 
system are polynomial functions of the e$, and requiring separation for all 
ti simultaneously sweeps out the full n{n + l)/2 space of symmetries and 
uniquely determines the nondegenerate potential. The infinitesimal distance 
in Jacobi elliptical coordinates Uj has the form 



4 trnZ =1 (^-e fc ) Ml= -4^~P(^ 



where -P(A) = n^ =1 (A — e^). However, it is well known that (1151) is a flat 
space metric for any polynomial P(X) of order < n and that each choice of 
such a P(X) defines an elliptic type multiplicative separable solution of the 
Laplace - Beltrami eigenvalue problem (with constant potential) in complex 
Euclidean n-space, |22j. The distinct cases are labeled by the degree of the 
polynomial and the multiplicities of its distinct roots. If for each distinct 
case we determine the most general potential that admits separation for all 
ti compatible with the multiplicity structure of the roots, we obtain a unique 
superintegrable system with nondegenerate potential and n{n + l)/2 second 
order symmetries, [2U I2"T] . These are the generic superintegrable systems. 
(Thus, for n = 3 there are 7 distinct cases for — ~ -P(A): 

(A- ei )(A-e 2 )(A-e 3 ), (A - e x )(A - e 2 ) 2 , (A - ei ) 3 , 

(A- ei )(A-e 2 ), (A- ei ) 2 , (A- ei ), 1, 

where e, 7^ Cj for i 7^ j. The first case corresponds to Jacobi elliptic coor- 
dinates.) The number of distinct generic superintegrable systems for each 
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integer n > 2 is ^2™ =0 p{j), where p(j) is the number of integer partitions of 
J- 

All of the generic separable systems, their potentials and their defining 
symmetries can be obtained from the basic Jacobi elliptic system in n dimen- 
sions by a complicated but well defined set of limit processes [2TJ EU [26] . In 
addition to these generic superintegrable systems there is an undetermined 
number of nongeneric systems. For n — 2 all the systems have been found, 
and now we give the results for n = 3. 

We review some of the details from reference [21] to show how each of the 
generic separable systems in three dimensions uniquely determines a nonde- 
generate superintegrable system that contains it. We begin by summariz- 
ing the full list of orthogonal separable systems in complex Euclidean space 
and the associated symmetries. (All of these systems have been classified, 
[22], and all can be obtained from the ultimate generic Jacobi elliptic co- 
ordinates by limiting processes [23, ETj.) Here, a "natural" basis for first 
order symmetries (Killing vectors) is given by p x = p x , p 2 = p y ,P3 = p z , 
■h — VPz — zp yi J 2 = zp x — xp z , J 3 = xp y — yp x in the classical case and 
Pi = d x , p 2 = d y , p 3 = d z , Ji = yd z - zd y , J 2 = zd x - xd z , J 3 = xd y - yd x 
in the quantum case. (In the operator characterizations for the quantum 
case, the classical product of two constants of the motion is replaced by the 
symmetrized product of the corresponding operator symmetries.) The free 
Hamiltonian is H Q = p\ + p\ + p\- In each case below we list the coordi- 
nates. The constants of the motion that characterize these coordinates can 
be found in [21] . We use the bracket notation of Bocher [26J to characterize 
each separable system. 



^ 2 _ j {u - e{){v - e-^jw - e{) 2 _ 2 (u - e 2 )(v - e 2 )(w - e 2 ) 



(&l - e 2 )(ei - e 3 ) ' (e 2 - e x ){e 2 - e 3 ) 

2 {u - e 3 )(v - e 3 )(w - e 3 ) 



z = c 



(e 3 - ei)(e 3 - e 2 ) 



[221] x 2 + y 2 = -c 2 



(u - ei)(v - e 1 )(w - ex) 



lei - e 2 



ei - e 2 



[(u - ei)(v - ei) + {u- e x ){w - e x ) + (v - e x )(w - e x )} 



(x - iy) 



2 = c 2 



(u - e x ){v - e x ){w - e x ) 



2 2 
Z = C 



(u - e 2 )(v - e 2 ){w - e 2 ) 



x — iy = -c 
y 2 



ei-e 2 (e 2 -ei) 2 

1 „ u 2 + v 2 + w 2 1 u 2 v 2 + u 2 w 2 + v 2 w 2 
uvw 2 u 3 v 3 w 3 
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1 uv uw vw 

z = —c( 1 1 ), x + iy = cuvw. 

2 w v u 

[311] x = Uu 2 + v 2 + w 2 + — + — + — ) + lc, 
4 u v w 2 

c (u 2 - l)(v 2 - l)(w 2 - 1) c (u 2 + l)(v 2 + l)(w 2 + 1) 

v = — A , Z = % - A • 

4 uvw 4 uvw 

UV UW VW 1/ 2 2 2\ 

32 x + zy = uvw, x — ly = — ( 1 1 ), z = -{u +v +w). 

w v u 2 

[41] x+iy = u 2 v 2 +v 2 w 2 +v 2 w 2 —-(u i +v 4 +w 4 ), x—iy = c 2 (u 2 +v 2 +w 2 ), z = licuvw. 
[5] x + iy = c(u + v + w), x — iy = -(u — v — w){u + v — w)(u + w — v) , 

z = —^(u 2 + v 2 + w 2 — 2{uv + uw + vw)). 

We summarize the remaining degenerate separable coordinates: 

Cylindrical type coordinates. All of these have one symmetry in 
common: L\ = p 2 . The 7 systems are, polar, Cartesian, light cone, elliptic, 
parabolic, hyperbolic and semihyperbolic. 

Complex sphere coordinates. These all have the symmetry L\ = 
Jf+Ji+Ji i n common. The 5 systems are spherical, horospherical, elliptical, 
hyperbolic, and semi-circular parabolic. 

Rotational types of coordinates. There are 3 of these systems, each of 
which is characterized by the fact that the momentum terms in one defining 
symmetry form a perfect square whereas the other two are not squares. 

In addition to these orthogonal coordinates, there is a class of nonorthog- 
onal heat-type separable coordinates that are related to the embedding of the 
heat equation in two dimensions into three dimensional complex Euclidean 
space. (22]- These coordinates are not present in real Euclidean space, only 
in real Minkowski spaces. The coordinates do not have any bearing on our 
further analysis as they do not occur in nondegenerate systems in three di- 
mensions. This is because they are characterized by an element of the Lie 
algebra pi + ip2 (not squared, i.e., a Killing vector) so they cannot occur for 
a nondegenerate system. 

Note that the first 7 separable systems are "generic," i.e., they occur in 
one, two or three - parameter families, whereas the remaining systems are 
special limiting cases of the generic ones. Each of the 7 generic Euclidean sep- 
arable systems depends on a scaling parameter c and up to three parameters 
e i> e 2, For each such set of coordinates there is exactly one nondegenerate 
superintegrable system that admits separation in these coordinates simulta- 
neously for all values of the parameters c,ej. Consider the system [23], for 
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example. If a nondegenerate superintegrable system separates in these co- 
ordinates for all values of the parameter c, then the space of second order 
symmetries must contain the 5 symmetries 

H = vl + p\ + p\ + V, S 1 = J 2 + J 2 + Jf + /i, S 2 = J 3 (Ji + U 2 ) + / 2 , 

#3 = (Px + «Pj,) 2 + /3, S4 = Pzfe + ^Py) + /4- 

It is straightforward to check that the 12 x 5 matrix of coefficients of the 
second derivative terms in the 12 Bertrand-Darboux equations associated 
with symmetries Si, • • • , S4 has rank 5 in general. Thus, there is at most one 
nondegenerate superintegrable system admitting these symmetries. Solving 
the Bertrand-Darboux equations for the potential we find the unique solution 

*r r / \ , 2 2 2n P 1 z 5(x 2 + y 2 -3z 2 ) 

V(x):= a (x +y +^) + ¥TW + ¥TW + {x + wY • 

Finally, we can use the symmetry conditions for this potential to obtain 
the full 6-dimensional space of second order symmetries. This is the su- 
perintegrable system III on the following table. The other six cases yield 
corresponding results. 

Theorem 5 Each of the 7 "generic" Euclidean separable systems determines 
a unique nondegenerate superintegrable system that permits separation simul- 
taneously for all values of the scaling parameter c and any other defining 
parameters Cj. For each of these systems there is a basis of 5 (strongly) func- 
tionally independent and 6 linearly independent second order symmetries. 
The corresponding nondegenerate potentials and basis of symmetries are: 

I [2111] V =^ + ^ + ^ + 5(x 2 + y 2 + z 2 ), (16) 
x 2 y 2 z 2 

Vi=pl+5x 2 + J i3 = (x iPxj - x jPxi ) 2 + a 2 -§ + a 2 -| , i>j. 

X i Xj Xj 

II [221] V = a(x 2 + y 2 + z 2 )+(3 7 ^- 3 + r ^— + - 2 , (17) 

(x + iyY (x + %y) 2 z 2 

S 1 = J-J + f 1 , S 2 =p 2 z + f 2l S 3 = J 2 3 +f 3 , 

S 4 = (Px + ipy) 2 + U, L 5 = ( J 2 - iJ x ) 2 + f 5 . 

2 , 2 , 2n , P l z , S(x 2 + y 2 -3z 2 ) 



III [23] V = a(x 2 + y 2 + z 2 ) + . ' . - + - '—— + 



x + iy) 2 (x + iy) 3 (x + iy) 4 

(1* 
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S 1 = J-J + f 1 , S 2 = ( J 2 - zJa) 2 + fa, S 3 = J 3 {J 2 - i JO + / 3 , 

5 4 = (P X + «Py) 2 + U, S 5 = + ZPy) + / 5 . 

IV [311] y = a(4x 2 + y 2 + z 2 )+(3x+^ + -^ } (19) 

^1 = Pi + fl, S 2 =p 2 y + f 2 , S 3 = P Z J 2 + fa, 

2 



5 4 = J 3 + /4, 5 5 = Jl + /5- 

V[32] V = a{Ax 2 + y 2 + z 2 ) + (3x + 7 + f * ' ~ '*] , (20) 

Si=p 2 + /\, S 2 = Jl + f 2 , S 3 = (p z -ip y )(J 2 + iJ 3 ) + fa, 
S 4 = Pz^2 - J3 + /*, 5 5 = (p z - ipy) 2 + / 5 . 

VI [41] V = a (z 2 - 2(z - iyf + A(x 2 + y 2 ))+/3 (2(x + iy) - 3(x - iy) 2 )+-f(x-i 

(21) 

«Si = (p x - ipj 2 + /1, S 2 =p 2 z + fa, S 3 =p z (J 2 + iJ 1 ) + fa, 

<?4 = J 3 (Ps - iPy) - ^(Px + iPyf + fa, S 5 = (J 2 + iJO 2 + &P-Jl + fa- 

3 1 13 

VII [5] V = a(x+iy)+p(-(x+iy) 2 +-z)+-f((x+iy) 3 + — (x-iy)+-(x+iy)z) 

4 4 lo 4 

(22) 

+S(—(x + iyf + — {x 2 + y 2 + z 2 ) + -(x + iy) 2 z), 
lo lo 8 

Si = (Ji + 2J2) 2 + HJ\{Px + ip y ) - Ji{Px + ip y ) + ^{P 2 y ~ P 2 Z ) ~ i-hPz + fl, 

S 2 = J 2 p z - J 3 p y + i(J 3 p x - JlPz) - -zPvP* + /2> 5 3 = (Pa + ^Py) 2 + fa, 

«5>4 = J 3 p z + UiPy + iJ 2 p x + 2Jip x + -p 2 + f 3 , S 5 = p z (p x + ip y ) + fa. 

In [21] we proved what was far from obvious, the fact that no other non- 
degenerate superintegrable system separates for any special case of ellipsoidal 
coordinates, i.e., fixed parameter. 

Theorem 6 A 3D Euclidean nondegenerate superintegrable system admits 
separation in a special case of the generic coordinates [2111], [221], [23], 
[311], [32], [41] or [5], respectively, if and only if it is equivalent via a Eu- 
clidean transformation to system [I], [II], [III], [IV], [V], [VI] or [VII], re- 
spectively. 
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This does not settle the problem of classifying all 3D nondegenerate su- 
perintegrable systems in complex Euclidean space, for we have not excluded 
the possibility of such systems that separate only in degenerate separable 
coordinates. In fact we have already studied two such systems in |19j : 

[O] V(x, y, z) = ax + (3y + 72 + S(x 2 + y 2 + z 2 ). 

1 3 
[00] V(x,y,z) = ^(x 2 + y 2 + -z 2 )+(3x + iy +-. (23) 
2 4 z z 

4 Polynomial ideals 

In this section we introduce a very different way of studying and classify- 
ing superintegrable systems, through polynomial ideals. Here we confine our 
analysis to 3D Euclidean superintegrable systems with nondegenerate poten- 
tials. Thus we can set G = in the 18 fundamental equations for the deriva- 
tives dia^ k . Due to the linear conditions (jHJ) all of the functions , B^ , C 1 - 7 
can be expressed in terms of the 10 basic terms ([9]). Since the fundamental 
equations admit 6 linearly independent solutions a hk the integrability con- 
ditions dia^ k = dia\ k for these equations must be satisfied identically. As 
follows from [19] . these conditions plus the integrability conditions (fTTl) for 
the potential allow us to compute the 30 derivatives deD^ of the 10 basic 
terms (equations ( 1601) in what follows). Each is a quadratic polynomial in the 
10 terms. In addition there are 5 quadratic conditions remaining, equation 
(31) in [T3] with G = 0. 

These 5 polynomials determine an ideal E'. Already we see that the 
values of the 10 terms at a fixed regular point must uniquely determine 
a superintegrable system. However, choosing those values such that the 5 
conditions /( a )-/( e ) ; listed below, are satisfied will not guarantee the existence 
of a solution, because the conditions may be violated for values of (2, y, z) 
away from the chosen regular point. To test this we compute the derivatives 
<9jE' and obtain a single new condition, the square of the quadratic expression 
listed below. The polynomial 1^' extends the ideal. Let E D E' be the 
ideal generated by the 6 quadratic polynomials, I^ a \ ■ ■ ■ , J^: 

j(o) = _^22 j g23 _j_ £23^33 + £12^13 + ^23 ^22 _ ypyp _ ^23^33 
j(b) _ (^33^2 + £12^33 _ ^33^22 _ ^12 ^33 _ ^13,^33 + ^12 ^22 

-B 12 A 22 + A 13 B 23 -(A 12 ) 2 

j(c) = ^23^33 + ^12^33 + (£12)2 + ^22^33 _ (£ >33 )2 _ ^12^33 _ (j g2 3) 2 
j(d) _ _ B 12 A 2Z _ ^33^23 + ^13 ^33 + ^12 ^23 
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j(e) _ _^j23^23 _|_ £<33^23 _|_ ^22^33 _ ^33^33 _j_ ^12^12 

/(/) = A 13 C <33 + 2A 13 jB 23 + jB 22 jB 33_ (jB 33 ) 2 + A 33 A 22_ (A 33 ) 2 

+ 2A 12 B 22 + (A 12 ) 2 - 25 12 A 22 + (5 12 ) 2 + B 23 C 33 - (B 23 ) 2 - 3(A 23 ) 2 . 

It can be verified with the Grobner basis package of Maple that <9jS C S, so 
that the system is closed under differentiation! This leads us to a fundamental 
result. 

Theorem 7 Choose the 10-tuple ([Sj) at a regular point, such that the 6 poly- 
nomial identities (Efy are satisfied. Then there exists one and only one 
Euclidean superintegrable system with nondegenerate potential that takes on 
these values at a point. 

We see that all possible nondegenerate 3D Euclidean superintegrable systems 
are encoded into the 6 quadratic polynomial identities. These identities de- 
fine an algebraic variety that generically has dimension 6, though there are 
singular points, such as the origin (0, ■••,(]), where the dimension of the 
tangent space is greater. This result gives us the means to classify all super- 
integrable systems. 

An issue is that many different 10-tuples correspond to the same super- 
integrable system. How do we sort this out? The key is that the Euclidean 
group E(3, C ) acts as a transformation group on the variety and gives rise 
to a foliation. The action of the translation subgroup is determined by the 
derivatives <9fc_D tJ that we have already determined (and will list below). 
The action of the rotation subgroup on the D %: > can be determined from the 
behavior of the canonical equations ^ under rotations. The local action on 
a 10-tuple is then given by 6 Lie derivatives that are a basis for the Euclidean 
Lie algebra e(3, C ). For "most" 10-tuples Do on the 6 dimensional variety 
the action of the Euclidean group is locally transitive with isotropy subgroup 
only the identity element. Thus the group action on such points sweeps out 
a solution surface homeomorphic to the 6 parameter E(3, C ) itself. This 
occurs for the generic Jacobi elliptic system with potential 




At the other extreme the isotropy subgroup of the origin (0, ■ • ■ , 0) is £"(3, C) 
itself, i.e., the point is fixed under the group action. This corresponds to the 
isotropic oscillator with potential 

V = a(x 2 + y 2 + z 2 ) + (3x + -fy + 5z. 
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More generally, the isotropy subgroup at D will be H and the Euclidean 
group action will sweep out a solution surface homeomorphic to the homo- 
geneous space E(3,C)/H and define a unique superintegrable system. For 
example, the isotropy subalgebra formed by the translation and rotation 
generators {Pi, P 2 , P3, J\ + iJ 2 } determines a new superintegrable system [A] 
with potential 

V = a((x- iy) 3 + Q(x 2 + y 2 + z 2 ))+/3 ((x - iy) 2 + 2(x + iy))+^(x-iy)+5z. 

Indeed, each class of Stackel equivalent Euclidean superintegrable systems 
is associated with a unique isotropy subalgebra of e(3, C ), although not 
all subalgebras occur. (Indeed, there is no isotropy subalgebra conjugate 
to {Pi, P 2 , P 3 }.) One way to find all superintegrable systems would be to 
determine a list of all subalgebras of e(3, C ), defined up to conjugacy, and 
then for each subalgebra to determine if it occurs as an isotropy subalgebra. 
Then we would have to resolve the degeneracy problem in which more than 
one superintegrable system may correspond to a single isotropy subalgebra. 

To begin our analysis of the ideal S we first determine how the rotation 
subalgebra so(3, C ) acts on the 10 variables ([9]) and their derivatives and 
decompose the representation spaces into so(3, C ) - irreducible pieces. The 
A % \ B l i and C 1 - 7 are 10 variables that, under the action of rotations, split 
into two irreducible blocks of dimension 3 and 7. 



X+1 = A 33 + 3B 12 - 2A 22 + i(3A 12 + B 33 + B 22 ) (25) 

X = -V2{C 33 + 2A 13 + B 23 ) (26) 

X_! = -A 33 -3B 12 + 2A 22 + i(3A 12 + B 33 + B 22 ) (27) 

y +3 = A 22 + 2B 12 + i(B 22 - 2A 12 ) (28) 

Y +2 = Vq(A 13 - B 23 + 2iA 23 ) (29) 

Y+1 = y^(3A 22 -2B 12 -4A 33 + t(B 22 -2A 12 -4B 33 )^j (30) 

Y = ^={2C 33 - A 13 -3B 23 ^ (31) 

YLi = ^(2B l2 + AA 33 -3A 22 + i(B 22 -2A l2 -AB 33 )) (32) 

Y_ 2 = V6(A 13 - B 23 - 2iA 23 ) (33) 

y_ 3 = -A 22 - 2B 12 + i(B 22 - 2A 12 ) (34) 
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Quadratics in the variables can also be decomposed into irreducible blocks. 
There are 2 one-dimensional representations, 3 of dimension 5, 1 of dimension 
7, 2 of dimension 9 and 1 of dimension 13. 

Zq 1ci ^ = Xq — 2X_iX + i (35) 



7 { 



(16) 



u i \ — / 

y 2 - 2y_ 1 y +1 + 2y_ 2 y +2 - 2y_ 3 y +3 (36) 

4 5 2 a) = Xli (37) 

45° = v^XoXi! (38) 

4 5a) = ^(^o 2 + ^-i^ + i) (39) 

4? = ^ - ^y y ±2 + ^y Tl y ±3 (40) 

4? - ^ r ±1 -|r Tl r ±2 + A w±3 (41) 

4- = $g-*Y_ lY+1 + ^Y» (42) 

4^ = -^Tl^±3 "I TFE^il^il 7=-^0^±2 (43) 

V 15 v 3 

4? = ^X ±1 F - ^X F ±1 + ^X T1 F ±2 (44) 

4 C) = -^X Y + ^X^Y +1 + ^X +1 Y^ (45) 
There is one 7-dimensional representation with highest weight vector 

Z%\ = X Y +3 - ^=X +1 Y +2 , (46) 



two 9-dimensional representations with highest weight vectors 

4? = y^- 2 -4y + iy + s m 



y +4 

and one 13-dimensional representation 

^(13) _ V 2 



4? = X +l Y + , (48) 



47 = • (49) 
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A linear combination of representations of the same dimension is another 
representation and if we define 

Z m = 2Z^-5ZW + 5Z£ C \ form = -2,-l,0,+l,+2. (50) 
W = %Z { Q la) - hZ { Q lb \ (51) 

the algebraic variety defining the nondegenerate superintegrable systems is 
given by 

Z m = W = for m = -2,-1,0, +1, +2. (52) 

If J x , J y and J z are Lie derivatives corresponding to rotation about the x, y 
and z axes, we define 

J + = iJ x + J y , J_ = i J x — J y and J3 = iJ z . 

then 

J+fm = \Z{l-m)(l+m+ l)/ m +i (53) 

J-fm = y/(l + m)(l-m + l)f m -i 

Jsfm Tnfm 

where f m is taken as one of X m , Y m , Z m or Wq. 

Derivatives of the X m and Y m are quadratics in these variables. The 
derivatives of the X m are linear combinations of the quadratics from the 
representations of dimensions 1 and 5. In particular, 

diXj G {2Z% a) + 5ZW : m = 0, ±1, ±2} U {Z ( Q 1A) } . (54) 

Hence the quadratic identities (15*21 can be used to write these derivatives as a 
sum of terms each of degree at least 1 in the X m . This means that whenever 
all of the X m vanish at a point, their derivatives also vanish and hence the 
set {X-i, Xo, X+x} is a relative invariant. 

The derivatives of the Y m are linear combinations of the quadratics from 
the representations of dimensions 5 and 9. 

diYj G {2Z^ a) +5Z£ b) : -2 < m < +2} U {hZf n a) -24Z% b) : -4 < m < +4} . 

(55) 

Hence they can be written as a sum of terms each of degree at least 1 in the 

Ymi SO 

Y-2, Y-!, Y , Y + i, Y +2 , Y +3 } 

is a relative invariant set. Note that from the dimension of the spaces con- 
taining the derivatives of the X m and Y m , there must be at least 3 linear 
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relations among the derivatives of the X m and 7 among the derivatives of 
the Y m . 

In a similar way we can find we can find relative invariant sets of quadrat- 
ics carrying a representation of the Lie algebra so(3, C ). For example, the 
following are relative invariant sets. 

R x = {X^,X ,X +1 }, (56) 

R2 = {Y-3, Y-2, y_i, Y , Y+i, Y +2 , Y +3 }, 

R 3 = {4Z£ 5a )-15^:m = 0,±l,±2}U{4 1A) }, 

R 4 = {3Z£ 5a )-5Z£ 5fe ):m = 0,±l,±2}U{4 1A) }, 

R 5 = {8Zt ) -5Z^:m = 0,±l,±2}, 

R 6 = R 5 U{5Z^ a) +6Z^ b) :m = 0,±l,±2,±3,±4}. 

Recall that the known superintegrable nondegenerate potentials are 
V, = a (x> + y> + S) + £ + l + ± (57) 

[x + iy) 6 [x + iy) A z 1 
Vni = a(x +y + z ) + p{x + iy) +- — + 



V IV = a(4x 2 + y 2 + z 2 )+(3x+^ + -^, 



(x + iy) 3 (x + iy) A 

y^ z 2 



V v = ct{Az + x + y ) + pz + — t + 



V VI = a(4x 2 + 4y 2 + z 2 - 2(x - iy) 3 ) + (3(2x + 2iy - 3(x - iy) 2 ) + >y(x-iy) + 



(x + iy) 2 (x + iy) 3 ' 

5_ 

Vyu = oi(x + iy) + P(3(x + iy) 2 + z) + 7(16(2; + iy) 3 + x — iy + 12z(x + iy)) 
+ 5(5(2; + iy) A + x 2 + y 2 + z 2 + 6(2; + iy) 2 z), 
V = a(x 2 + y 2 + z 2 ) + f3x + ^y + Sz, 

Voo = a(Ax 2 + Ay 2 + z 2 )+ (3x + -fy+ —, 
V A = a{{x-iy) 3 + 6{x 2 + y 2 + z 2 ))+l3({x-iy) 2 + 2x + 2iy)+-f{x-iy) + 5z. 

The correspondence between relative invariant sets and potentials is in the 
accompanying table. 
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The action of the Euclidean translation generators on the 10 basis mono- 
mials can also be written in terms of the irreducible representations of so(3, C ). 
(Indeed these equations are much simpler than when written directly in terms 
of the A ij , B ij , C ij .) Using the notation 



«9 4 



id y ± d x 



Z(f ) = 5Z^ + 2Z% a) , ZfP = - 5Z^ . 

we obtain the fundamental differential relations: 



d-X +1 
d z X +1 
d + X 

d z X_ x 



1 



r(5X) 



(1A) 



30V6 



d + X +1 = -U 



30 



(5X) 
+2 > 



7 (5X) 



30^ 



1 7 (5X) 

30V2 +1 ' 



d z X = 



r(5X) 



r(lA) 



9V2 



30 

60 - 1 ' 



30^ 



zf x) - -z 



(1A) 



(58) 
(59) 



(60) 



d-Y +3 
d z Y +3 
d+Y +2 
d-Y +1 



_i 7 (9y) , 1 7 (5X) 

lSOv/7 +2 35 Z+2 ' 

r(9V) 



d + Y_ 



+3 



180^ + 



60V21 



90 

7 (9Y) 



J +4 • 



+ 



35^' 



7 (bX) 



1 Z (9Y) 

1 



30V42 



d z Y +2 = 
V2 



1 



Z 



(9Y) 
+2 



Z (9Y) Z (5X) 



35v/5 



30V42 



+ 



1 7 (5X) 

35^ +2 ' 
1 7 (sy) 



(61) 



1 



12V105 



35V15 



(5X) 
+2 5 
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d z Y +1 = ^i-Lzjf, d _ Yo = ^—2FP + ^-=X<™ , 

+ 12^210 +1 35^ +1 18^ _1 35v/5 1 



_i 7 (9Y) 1 y (5X) ~ v _ 1 7 (9Y) y (5X) 



12V/T05 ~ 2 35^ " 2 ' + ^ 30^42 35^ ° : 



d,Y. x = 1=2? + ^ 9_Y1 2 = ^Zi 9 3 y) , 

12V/2T0 35^ 60^ 

= -J= z^ + -4=z^, d,Y. 2 = — 1 z?P + ^z% x \ 

60^ 35^ 30^ 35^6 



3 90 ~ 4 ' + 3 180^ ~ 2 35 - 2 
d z Ys = -~^=Z^P. 



I8OV2' 



In the following table we describe each of the known superintegrable sys- 
tems in terms of variables adapted to the rotation group action, for this it is 
convenient to choose the 10 constrained variables in the form X i} i = 1 . . . 3 
and Yj, j = 1 . . .7 with dx and dy, respectively, as the number of indepen- 
dent variables on which these variables depend. These are defined by 

X 1 = 2A 13 + B 23 + C 33 = -4^, X 2 = 2A 22 - A 33 - 3B 12 = X ~ 1 ~ X+1 , 

v 2 2 

X 3 = 3A 12 + B 33 + B 22 = X - 1+ 2 X+ \ Yi = \{Y + , - F_ 3 ), 



(62) 



y 2 = i(r +3 + r-), f 3 = ^(r +2 - r_ 2 ), y 4 = ^=(r +2 + r_ 2 ), 



(63) 
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[Xi,X 2 , x 3 \ 



d x 



\Y U Y 2 ,Y 3 ,Y A ,Y 5 ,Y 6 ,Y 7 ] 



9_ , _9_ , _9_ 

x 2 y 2 z 2 



V, 



11 



V, 



III 



V, 



IV 



v 



VI 

Vvii 
Vo 

Voo 
Va 



_9_ _j_ _9_ 

y 2 ~r z 2 



SL 

,2 



3 _3 3 

x ' y ' 2 



6 __6i_ 3 

x+iy ' x+iy ' z 



9i 



x+iy" 1 x+iy 



0-3 3 

' 2/ ' z 



6 6i 



.0 



x+iy ' x+iy 

[0,0,0] 
[0,0,0] 

[0,0,0] 



,0 



3 
3 

2 
3 

1 

3 

2 
2 

1 
2 

1 
1 


1 




1 
1 






3 _3 q g _3 _3 _6 

x ' y ' ' ' x ' y ' z 



6(x-iy) 6i(x-iy) p, p, 6 
" (x+iy) 2 ' (x+iy) 2 ' U ' U ' 



6i _6 
x+iy' 1 x+iy' 1 z 



6(x 2 +y 2 -2z 2 ) 6i(x 2 +y 2 -2z 2 ) 6iz 
(x+iy) 3 ' (x+iy) 3 ' (x+iy) 2 ' 

62 6 6i g 

(x+iy) 2 ' x+iy' x+iy" 1 \ 

0,-^, o,o,o,--,--l 



6(x-iy) 6i(x~iy) n , , (j 
(x+iy) 2 ' (x+iy) 2 



-,0,0 



x+iy' x+iy 



.0 



[6,-6i,0,0,0,0,-f] 
[-48(x + iy), -48i(x + iy), 12i, 12, 0, 0, 0] 
[0,0,0,0,0,0,0] 
[0,0,0,0,0,0,-f] 
[-2,2^,0,0,0,0,0] 



In principle one could classify all possibilities by referring to distinct cases 
exhibited in the accompanying table. Here, however, we use the preceding 
algebraic and differential conditions, together with the coordinates in which 
the corresponding nondegenerate system could separate, to demonstrate that 
our 10 known superintegrable systems are the only ones possible. 
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5 Completion of the proof 



We know that in addition to the generic superintegrable systems, the only 
possible superintegrable systems are those that are multiseparable in non- 
generic coordinates. Our strategy is to consider each nongeneric separable 
system in a given standard form and use the integrability conditions associ- 
ated with the corresponding separable potential. If a superintegrable system 
permits separation in these coordinates, then by a suitable Euclidean trans- 
formation, we can assume the system permits separation in this standard 
form. This information is then used together with the six algebraic condi- 
tions l( a \ • • • , (124 j) . to deduce all the information available from algebraic 
conditions. At that point the differential equations (I6UI) for the D u can be 
solved in a straight forward manner to obtain the final possible superinte- 
grable systems. In some cases the algebraic conditions alone suffice and the 
differential equations are unnecessary. We proceed on a case by case basis. 

5.1 Cylindrical systems 

For cylindrical- type systems the potential splits off the z variable , i.e., the 
potential satisfies V13 = 0, V23 = in equations ([31). This implies that A 13 = 
B 13 = C 13 = and A 23 = B 23 = C 23 = 0. From the equations for X h (i = 
1, 2, 3) and Yj, (j = 1, • • • , 7) we can deduce that Y 7 = —2X 3 . In addition it 
is also easy to conclude that Y3 = I4 = and X\ = Y5, X 2 = Yg. 

If we add the requirement of Cartesian coordinate separation then A 12 = 
£12 = c i2 = Q jf x 3 = we obtain potential V . If X 3 ^ then X 3 = 3/z 
If Xi = X 2 = then we have potential Voo- If one of X±, X 2 is not zero this 
leads directly to potential Vj. 

For separation in cylindrical coordinates x = r cos 8, y = r sin 9, z, the 
following conditions must apply: 

v xz = 0, v yz = 0, 

(x 2 - y 2 )V xy + xy(Vyy - V xx ) + 3xV y - 3yV x = 0. 

The last condition is equivalent to do(rd r (r 2 V)) = where r 2 = x 2 + y 2 . 
Solving the algebraic conditions that result, we determine that 

X x = Y 5 = -G(l + 4) - -, X 2 = Y 6 = G{- + V -)-- : 

x A x y x y 

Y x = G(-3 + ^) + ~ Y 2 = G(- - 3 V -) - Y 3 = Y A = 0, 

x z x y x y 
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where G is an unknown function. In addition we deduce that Y 7 = — 2X 3 . 
It is then easy to show from the differential equations that X 3 = | or and 
that G — 0. We conclude that separation of this type occurs in cases Vj and 
V IV . 

For parabolic cylinder coordinates x — |(£ 2 — i] 2 ), y = £77, z, the condi- 
tions on the potential have the form 

V xz = 0, V yz = 2xV xy + y{V yy - V xx ) + W y = 0. 

This implies that 

X, = -2F, X 2 = 2-F - -, X 3 = -C, 

y y 

Y 1 = -2F, Y 2 = 2-F --,Y 3 = Y 4 = 0, 

y y 

Y 5 = -2F, Y 6 = 2-F --,Y 7 = 2C. 

y y 

The remaining differential equations require that F = and C — |. This 
type occurs in case Vjy. 

For elliptic cylinder coordinates x = cosh A cos B,y = sinh A sin B, z, the 
integrability conditions for the potential have the form 

V zx = 0, V vz = 0, (x 2 -y 2 - l)V xy + xy(V yy - V xx ) + 3(xV y - yV x ) = 0. 

This and the algebraic conditions imply 

X 1 = (- + - + -)G--, X 2 = (-l-^ + \)G--, X 3 = -C, 

y x xy x y z y z y 

Y 1 = (3--y--—)G + -, Y 2 = (~ + 3 + \)G--, Y 3 = Y A = 0, 
y x xy x y z y z y 

Y 5 ^(- + y + -)G--, Y 6 = (-l-^ + - 2 )G--, Y 7 = 2C. 

y x xy x y z y A y 

The remaining differential equations require G — 0, and C = — | or corre- 
sponding to systems Vi and Viy. 

In semihyperbolic coordinates x + iy = 4i(u + v), x — iy = 2i(u — v) 2 the 
extra integrability condition is 

(1 + ix + y)(V xx - V yy ) + 2(-2i - x + iy)Ks/ + 3iV x - 3V y = 0. 
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The algebraic conditions yield the requirements 

Xi — Y 5 — G, X 2 = —G, X 3 = — C, Y 3 = Y 4 = 0, 

3 3 1 

Y x = -i+ % -{x- iy)G, Y 2 = -- + -(-x + iy)G, Y 6 = iG, Y 7 = 2C. 

This leads to potentials Va and Vyj. 

For hyperbolic coordinates x + iy = rs,x — iy = (r 2 + s 2 )/rs,z, the 
integrability condition is 

(1 + ixy){V yy - V xx ) + i[x 2 -y 2 - 2)V xy + 3i(xV y - yV x ) = 0. 

The algebraic conditions imply Y 7 = 2X 3 = 2C and 

X 1 = Y 5 = ( xy -iy 2 -2i)G--^—, X 2 = Y 6 = -( X 2 -t X y-2)G--^-, 

x + iy x + iy 

3yx 2 — 2ix — y 3 — 2y 6(x — iy) x 3 — 3xy 2 — 2x + 2iy ,6(x — 

1 = r~ — 1 • \2 ' 2 = r~ g—%- — - 

x + iy [x + iy) A x + iy [x + 

This yields potential V/j. 

5.2 Radial- type coordinates 

We consider systems that have a radial coordinate r as one of the separable 
coordinates. The two other coordinates are separable on the complex two 
dimensional sphere. We first consider spherical coordinates x = r sin 9 cosy?, 
y = rsin^siny?, z = rcosO. The integrability conditions on the potential 
have the form 

(x 2 - y 2 )V xy + xzV yz - yzV xz + xy(V yy - V xx ) + 3xV y - 3yV x = 0, 

{x 2 - z 2 )V xz + xz(V zz - V xx ) + xyV yz - zyV xy + 3xV z - 3zV x = 0, 
{y 2 - z 2 )V yz + yz(V zz - V yy ) + xyV xz - zxV xy + 3yV z - 3zV y = 0, 

xV yz - yV xz = 0. 

Note that the first three conditions are not independent and only two are 
required. For any potential that separates in spherical coordinates, one ad- 
ditional condition is required. Indeed, if r, u and v are any form of separable 
spherical-type coordinates then the potential must have the functional form 

V = f(r)+g(u,v)/r 2 , (64) 
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it being understood that u and v are coordinates on the complex two dimen- 
sional sphere and r is the radius. It is then clear that r 2 V = r 2 f(r) + g(u, v). 
As a consequence there are the conditions d r d\(r 2 V) = , where X = u,v. 
Noting that 

xd x F + yd y F + zd z F = DF = rd r F 

and that 

J\F = yd z F — zd y F = a(u, v)d u F + b(u, v)d v F, 

with similar expressions for J 2 F and J3F, we conclude that the conditions 
( Id^|) are equivalent to any two of the three conditions ^JiD{r 2 V) = 0. These 
are indeed the three conditions we have given. If we now solve all the alge- 
braic conditions, we determine that 



(x 2 + y 2 ) 3 _ (x 2 + y 2 ) 3 3 6 

A x — r 5 — Lt , A 2 — t, <jt , ^3 — _ , r 7 — — i 

xy x y z y z z 

xy x y z y 

From this we see that the remaining differential equations give G = and 
we obtain solution Vj. 

We now consider horospherical coordinates on a complex 2 sphere, viz. 

• r / 2 2x . r . u 

x + ly = —i—lu + v ), x — ly = 1—, z = —ir — . 

V V V 

The extra integrability condition in this case is 

z(V xx - V yy ) + 2izV xy - (x + iy){V xz + iV yz ) = 0. 
Solving the algebraic conditions we conclude that 

x x = ix 2 = (^±m g - x 3 = {x+ ; y)2 G+ 3 -, 

z x + iy z 2 z 

Yx = iY 2 = -4 . - , G - &~ lV } 2 , Y 3 = iY, = -2iG, 
[x + iy) [x + iy) A 

Y 5 = i Y 6 = -A^l G - r ^, Y 7 = -2^±^G- 6 -. 

z [x + iy) z l z 

The derivative conditions give G = 0, so this corresponds to solution V/j. 
Conical coordinates are also radial-type: 

2 2 (tt-e x )(v-ei) 2 2 (u-e 2 )(v-e 2 ) 
x = r — -, y = r 



(ei - e 2 ) (ei - e 3 ) (e 2 - e 1 )(e 2 - e 3 ) 
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2 2 (u-e 3 )(u-e 3 ) 

,0 — / 



(e 3 - e 2 )(e 3 - d) 
The extra integrability condition is 

3(e 2 -e 3 )^K+3(e 3 -ei)x2;V^+3(ei-e2)^V;+^[(e2-e3)K a; +(e 3 -ei)V^j / +(ei-e2)V; 2 

+z[(e 3 -e 1 )y 2 +(e 2 -e 3 )x 2 +(e2-e 1 )z 2 ]V X y+y[(e 1 -e 2 )z 2 +(e2-e 3 )x 2 +(e 1 -e 3 )y 2 ]V xz 

+x[(e 1 - e 2 )2; 2 + (e 3 - e 2 )x 2 + (ei - e 3 )y 2 ]V^ = 0. 
The algebraic conditions yield immediately solution Vj with 

x - 3 x - 3 x - 3 y- 3 y- 3 

Ai — , A2 — , A 3 — -, Ii — — , 12 — , 

x y z x y 

Y 3 = Y A = 0, Y 5 = --, Y 6 = --, F 6 ? = --. 

x y x 

For degenerate type elliptic polar coordinates (type 1) we can write 

r cosh A sinh£> 

x + zy = — — -, 2x = r\ ; — r\, z = r tanh Atanh±5. 

cosh A cosh B cosh is smh A 

The extra integrability condition is 

3(x+iy) 2 V z -3xzV x -3i(2x+iy)zV y -2i(x+iy)(z 2 +ixy)V yz -2(y 2 +z 2 )(x+iy)V xz 

+2iz(z 2 + y 2 )V xy + z(x + iy) 2 V zz + z(z 2 + y 2 )^ - z(x 2 + z 2 + 2ixy)V yy = 0. 
Solving the algebraic conditions we deduce that 

2 6 2 6i 

Xi = — (y-ix)G — , X 2 = -i-(y 2 -x 2 + z 2 -ixy)(y-ix)G — , 

x x + ly x x + ty 

X 3 = - — {z 2 +y 2 -txy){y-tx) 2 G+-, Y 1 = --(-y 3 +3x 2 y+2z 2 y-6tz 2 x)G-6 ^ ~ %V > 



XZ Z X [X + iy) 

Y 2 = --(-3ixy 2 +ix 3 +2z 2 y-Qiz 2 x)G-Qi . - - W , , Y 3 = iY A = 2z <V ~ G, 
x (x + iy) 1 x 

3 6 
F 5 = — {-3y 2 + 5ixy + 2z 2 ){y -ix)G - 
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x ' x + iy 

Y 6 = — (-8y 2 +13ixy+3x 2 +2z 2 )(y-ix)G , Y 7 = (-2y 2 +2ixy+3z 2 )(y-ix) 2 G — . 

6 x + iy xz z 

The differential conditions require G = 0, leading to a type Vu potential. 
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For degenerate elliptic coordinates (type 2) on the complex 2 sphere we 
have 

/ 2 2\2 3 3 % 11 

x + %y = ruv, x — %y = -r(u + v ) u v , z — — r . 

4 2 uv 

The corresponding integrability condition is 

3(z 2 + ixy — y 2 )V x + 3i{z 2 — x 2 — ixy)V y — 3iz(y — ix)V z 

—i{—ixy 2 + y 3 + iz 2 x + yz 2 )V xx + i(ix 3 — x 2 y + iz 2 x + tz 2 )V xx 
+i(—ix + y)(x 2 + y 2 )V zz + 2(x 2 y + yz 2 + ixy 2 + ixz 2 )V xy 
-2iz{x 2 + y 2 )V yz - 2z{x 2 + y 2 )V xz = 0. 
The solutions to the algebraic conditions are 

9 9i 
Xi = -2iz(ix + 2y)(y-ix)G -, X 2 = 2z{y + ix)(y - ix)G - 



x + iy ' x + iy' 

Qiz 



X 3 = 2(-ix + y){x 2 + y )G, Y 3 = 2i(yz 2 + iz 2 x - ixy 2 - x 2 y)G + 



(x + iy) 2 



. (-3y 2 - 3x 2 + Az 2 )z{ix + y) {-x 2 - y 2 + 2z 2 ) 

Yi = -1^2 = ^ : G + 6 — — — , 

ix — y [x + iy) 6 

6 6 
F 4 = {2yz 2 +2iz 2 x— y 3 +ixy 2 +x 2 y— ix 3 )G+- — -, Y 5 = iz(y+3ix)(y—ix)G+ 



(x + iy) 2 ' x + iy' 



6i 

Y 6 = -z{ix + 3y)(y - ix)G + ——, Y 7 = (-ix + y)(x 2 + y 2 )G. 

x + iy 

The differential conditions hold only if G = 0. This is system V in . 



5.3 Spheroidal coordinates 

We take these as 

x = sinh A cos B cos <p, y — sinh A cos B sin ip, z = cosh A sin B. 
The integrability conditions for the potential are 

-3zV x + 3xV z + zx{V zz - V xx ) - zyV xy + (1 + x 2 + y 2 - z 2 )V zx = 0, 

-3zV y + 3yV z + zy(V zz - V yy ) - zxV xy + (1 + x 2 + y 2 - z 2 )V zy = 0, 

yV zx - xV zy = 0. 
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The solutions of the algebraic conditions are 

X 1 = Y 5 = -V(x 2 + y 2 )G--, X 2 = Y 6 =(x 2 + y 2 )G--, X 3 = -, 

xx y z 

Y, = - V -{-y 2 + 3x 2 )G + -, Y 2 = (-3y 2 + x 2 )G --,Y 7 = --. 
x rx y z 

From the differential conditions we see that G = 0, and obtain potential Vj. 



5.4 Horospherical coordinates 

These are 

x + iy = v /pz7, x-iy = 4 — , z = 1yfpv\L. 



The corresponding integrability conditions for the potential are 

(x 2 -ixy-z 2 )V zx +(yx-iy 2 +iz 2 )V zy +i(x+iy)zV xy +zx(V zz -V x ^ = 0, 

(x 2 - y 2 )V xy + xy{V yy - V xx ) + zxV zy - yzV zx - 3yV x + 3xV y = 0, 
z{V xx - V yy ) - 2izV xy + (ix + y)V zy + (-x + iy)V zx = 0. 
The solutions to all the algebraic conditions are 



i(x + iy) 6 i(x + iy) 2 3 
Ai = -%X 2 = G — - , A 3 = - G + -, 

z x + ly z z z 



Y-v = -iY 2 = G _ 6-^-5-, Y 3 = iY 4 = 2G, 

x + ly [x + lyy 

Y 5 = -,r 6 = A i{x + iy) G - Y 7 = -2i^±pl G - 6 -. 

z x + iy z 2 z 

The differential conditions require G = and this gives potential V/j. 
5.5 Rotational parabolic coordinates 

For these coordinates x = £77 cos ip , y = £77 simp, z = |(£ 2 — rj 2 ). The 
required conditions on the potential are 

xy(Vyy - V xx ) + (x 2 - y 2 )V xy - yzV yz + xzV^ - 3yV x + 3x1^ = 0, 

x 2 {V xx - V zz ) + y 2 {V yy - V zz ) + 2xyK y + 2zxV yz + 2xzV zx + 3xV x + 3yV y = 0, 

xV zy - yV zx = 0. 
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These integrability conditions directly produce the solution 

X% = , X2 = — , X3 = 0, Y\ — — , Y2 = — , 

x y x y 

Y 3 = Y, = 0, Y 5 = --, Y 6 = --, Y 7 = 0. 

x y 

This is a permuted version of potential Viy. 

We have covered all possibilities for separable coordinates and found ex- 
actly which superintegrable system separates in each coordinate system It 
follows that our list of 10 superintegrable systems is complete. Another in- 
teresting consequence of this analysis is 

Theorem 8 For every orthogonal separable coordinate system there is at 
least one nondegenerate superintegrable system that separates in these coor- 
dinates. 

On the other hand, no nondegenerate superintegrable system permits sep- 
aration in nonorthogonal heat-type coordinates. Potential Vvn is the only 
generic system that separates in generic coordinates alone. 

6 Outlook 

The basic structure and classification problems for 2D second order super- 
integrable systems have been solved, [El EHl EH [30, l3l] . For 3D systems 
the corresponding problems are much more complicated, but we have now 
achieved a verifiably complete classification of the possible nondegenerate 
potentials in 3D Euclidean space. There are 10 such potentials, as compared 
to 11 in 2D. To finish the classification of nondegenerate potentials for all 3D 
conformally flat spaces the main task remaining is the classification on the 
3-sphere, probably not difficult. This is because all conformally flat systems 
can be obtained from flat space and the 3-sphere by Stackel transforms. The 
new idea used here that made the complete verifiable classification practical 
was the association of nondegenerate superintegrable systems with points on 
an algebraic variety on which the Euclidean group acts to produce foliations. 
In the future we hope to refine this approach to give a direct classification 
using only the algebraic variety and group action. Here we had also to rely 
on basic results from separation of variables theory to simplify the calcula- 
tions. In distinction to the 2D case, which is special, the 3D classification 
problem seems to have all of the ingredients that go into the corresponding 
nondegenerate potential classification problem in n dimensions. The number 
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of nondegenerate potentials grows rapidly with dimension: the number of 
generic potentials alone is ^™ =0 p(j), where p(j) is the number of partitions 
of j. The algebraic variety approach should be generalizable to this case. 

Nondegenerate potentials for 3D superintegrable systems are just the 
most symmetric. There is also "fine structure," a hierarchy of various classes 
of degenerate potentials with fewer than 4 parameters. The structure and 
classification theory for these systems has just gotten underway, with initial 
results for 3 parameter FLI systems. [20] ■ Sometimes a quadratic algebra 
structure exists and sometimes it does not. Extension of these methods 
to complete the fine structure analysis for 3D systems appears relatively 
straightforward. The analysis can be extended to 2 parameter and 1 param- 
eter potentials with 5 functionally linearly independent second order symme- 
tries. Here first order PDEs for the potential appear as well as second order, 
and Killing vectors may occur. Another class of 3D superintegrable systems 
is that for which the 5 functionally independent symmetries are functionally 
linearly dependent. This class is related to the Calogero potential [32l I33~t 134"] 
and necessarily leads to first order PDEs for the potential, as well as second 
order |9] . However, the integrability methods discussed here should be able to 
handle this class with no special difficulties. On a deeper level, we think that 
the algebraic geometry approach can be extended to determine the possible 
superintegrable systems in all these cases. 

Finally, the algebraic geometry related results that we have described in 
this paper suggest strongly that there is an underlying geometric structure 
to superintegrable systems that is not apparent from the usual presentations 
of these systems. 

Acknowledgement: The authors wish to thank Thomas Wolf and Greg 
Reid for very helpful consultations on Grobner basis techniques and on nu- 
merical methods in the study of algebraic varieties. 

References 

[1] Wojciechowski S., Superintegrability of the Calogero- Moser System. 
Phys. Lett, V. A 95, 279-281, (1983). 

[2] Evans N.W., Superintegrability in Classical Mechanics; Phys. Rev. V. 
A 41, 5666-5676 (1990); Group Theory of the Smorodinsky-Winternitz 
System; J. Math. Phys. V. 32, 3369, (1991) 

[3] Evans N.W., Super-Integrability of the Winternitz System; Phys. Lett. 
V.A 147, 483-486, (1990). 



32 



[4] Fris J., Mandrosov V., Smorodinsky Ya.A, Uhlir M. and Winternitz 
P., On Higher Symmetries in Quantum Mechanics; Phys. Lett. V.16, 
354-356, (1965). 

[5] Fris J., Smorodinskii Ya.A., Uhlir M. and Winternitz P., Symmetry 
Groups in Classical and Quantum Mechanics; Sov. J. Nucl. Phys. V.4, 
444-450, (1967). 

[6] Makarov A. A., Smorodinsky Ya.A., Valiev Kh. and Winternitz P., A 
Systematic Search for Nonrelativistic Systems with Dynamical Symme- 
tries. Nuovo Cimento, V. 52, 1061-1084, (1967). 

[7] Calogero F., Solution of a Three-Body Problem in One Dimension. J. 
Math. Phys. V.10, 2191-2196, (1969). 

[8] Cisneros A. and Mcintosh H.V., Symmetry of the Two-Dimensional Hy- 
drogen Atom. J. Math. Phys. V.10, 277-286, (1969). 

[9] Kalnins E.G., Kress J.M, and Miller W.Jr., Second order superintegrable 
systems in conformally flat spaces. V: 2D and 3D quantum systems. J. 
Math. Phys., V.47, 093501, (2006). 

[10] L. G. Mardoyan, G. S. Pogosyan, A. N. Sissakian and V. M. Ter- 
Antonyan. Elliptic Basis for a Circular Oscillator. Nuovo Cimento, B 
88, 43 (1985). Two-Dimensional Hydrogen Atom: I. Elliptic Bases; 
Theor. Math. Phys. 61, 1021 (1984); Hidden symmetry, Separation of 
Variables and Interbasis Expansions in the Two-Dimensional Hydrogen 
Atom. J. Phys., A 18, 455 (1985). 

[11] Ya. A. Granovsky, A. S. Zhedanov and I. M. Lutzenko. Quadratic Al- 
gebra as a 'Hidden' Symmetry of the Hartmann Potential; J. Phys. A 
24, 3887 (1991). 

[12] Bonatos D., Daskaloyannis C. and Kokkotas K., Deformed Oscillator Al- 
gebras for Two-Dimensional Quantum Superintegrable Systems; Phys. 
Rev., V.A 50, 3700-3709, (1994). 

[13] Grosche C, Pogosyan G.S., Sissakian A.N., Path Integral Discussion for 
Smorodinsky - Winternitz Potentials:!. Two - and Three Dimensional 
Euclidean Space. Fortschritte der Physik, V. 43, 453-521, (1995). 

[14] Kalnins E.G., Kress J.M., Miller W.Jr. and Pogosyan G.S., Complete- 
ness of superintegrability in two-dimensional constant curvature spaces. 
J. Phys. A: Math Gen., V.34, 4705-4720, (2001). 



33 



[15] Kalnins E.G., Kress J.MN., and Winternitz P., Superintegrability in a 
two-dimensional space of non- constant curvature. J. Math. Phys. V.43, 
970-983, (2002). 

[16] Kalnins E.G., Kress J.M., Miller, W.Jr. and Winternitz P., Superin- 
tegrable systems in Darboux spaces. J. Math. Phys., V.44, 5811-5848, 
(2003). 

[17] Ranada M.F., Superintegrable n=2 systems, quadratic constants of mo- 
tion, and potentials of Drach. J. Math. Phys., V.38, 4165-78, (1997). 

[18] Kalnins E.G., Miller W.Jr., Williams G.C. and Pogosyan G.S., On su- 
perintegrable symmetry-breaking potentials in n-dimensional Euclidean 
space. J. Phys. A: Math Gen., V.35, 4655-4720, (2002). 

[19] Kalnins E.G., Kress J.M, and Miller W.Jr., Second order superintegrable 
systems in conformally flat spaces. Ill: 3D classical structure theory. J. 
Math. Phys., V.46, 103507, (2005). 

[20] Kalnins E.G., Kress J.M, and Miller W.Jr., Fine structure for second or- 
der superintegrable systems (to appear in the procedings volumes of the 
IMA program on " Symmetries and Overdetermined Systems of Partial 
Differential Equations." ) 2007. 

[21] Kalnins E.G., Kress J.M, and Miller W.Jr., Second order superinte- 
grable systems in conformally flat spaces. IV: The classical 3D Stackel 
transform and 3D classification theory. J. Math. Phys., V.47, 043514, 
(2006). 

[22] Kalnins E.G., Separation of Variables for Riemannian Spaces of Con- 
stant Curvature, Pitman, Monographs and Surveys in Pure and Applied 
Mathematics V.28, 184-208, Longman, Essex, England, 1986. 

[23] Miller W.Jr., Mechanisms for variable separation in partial differential 
equations and their relationship to group theory. In Symmetries and 
Non-linear Phenomena pp. 188-221, World Scientific, 1988. 

[24] E.G. Kalnins, J.M.Kress, W.Miller Jr. and G.S. Pogosyan. Non degener- 
ate superintegrable systems in n dimensional complex Euclidean spaces. 
(To appear in Physics of Atomic Nuclei) 

[25] Kalnins E.G., Kress J.M, and Miller W.Jr., Second order superintegrable 
systems in conformally flat spaces. I: 2D classical structure theory. J. 
Math. Phys., V.46, 053509 (2005). 



34 



[26] M. Bocher. Uber die Riehenentwickelungen der Potentialtheory. Sym- 
metry and Separation of Variables. Teubner, Leipzig, 1894. 

[27] E. G. Kalnins, W. Miller and G. K. Reid, Separation of variables for 
Riemannian spaces of constant curvature. I. Orthogonal separable coor- 
dinates for S c and E nC . Proc R. Soc. Lond. A 39, 183-206 (1984). 

[28] Kalnins E.G., Kress J.M, and Miller W.Jr., Second order superinte- 
grable systems in conformally flat spaces. II: The classical 2D Stackel 
transform. J. Math. Phys., V.46, 053510, (2005). 

[29] E. G. Kalnins, W. Miller, Jr. and G. S. Pogosyan. Completeness of 
multiseparable superintegrability in E 2t c- J- Phys. A: Math Gen. 33, 
4105 (2000). 

[30] E. G. Kalnins, W. Miller Jr. and G. S. Pogosyan. Complete- 
ness of multiseparable superintegrability on the complex 2-sphere. 
J. Phys. A: Math Gen. 33, 6791-6806 (2000). 

[31] C. Daskaloyannis and K Ypsilantis. Unified treatment and classification 
of superintegrable systems with integrals quadratic in momenta on a two 
dimensional manifold. J. Math. Phys. 47, 042904 (2006). 

[32] F. Calogero. Solution to the one-dimensional iV-body problems with 
quadratic and/or inversely quadratic pair potentials. J. Math. Phys. 12, 
419-436, (1971). 

[33] S. Rauch-Wojciechowski and C. Waksjo. What an effective criterion of 
separability says about the Calogero type systems. J. Nonlinear Math. 
Phys. 12, Suppl. 1 535-547, (2005). 

[34] J. T. Horwood, R. G. McLenaghan and R. G. Smirnov. Invariant clas- 
sification of orthogonally separable Hamiltonian systems in Euclidean 
space. Comm. Math. Phys. 259, 679-709, (2005); 



35 



